Non-Markovian master equation for a damped oscillator with time-varying parameters 
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We derive an exact non-Markovian master equation that generalizes the previous work [Hu, Paz 
and Zhang, Phys. Rev. D 45, 2843 (1992)] to damped harmonic oscillators with time-varying 
parameters. This is achieved by exploiting the linearity of the system and operator solution in 
Heisenberg picture. Our equation governs the non-Markovian quantum dynamics when the system 
is modulated by external devices. As an application, we apply our equation to parity kick decoupling 
problems. The time-dependent dissipative coefficients in the master equation are shown to be 
modified drastically when the system is driven by it pulses. For coherence protection to be effective, 
our numerical results indicate that kicking period should be shorter than memory time of the bath. 
The effects of using soft pulses in an ohmic bath are also discussed. 

PACS numbers: 03.65.Yz, 42.50.Lc 
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I. INTRODUCTION 

An oscillator linearly coupled with a harmonic oscilla- 
tor bath has been an important model for studying quan- 
tum dissipation and decoherence. While the solution 
for the damped harmonic oscillator is well-known under 
Born- and/or Markov-approximation [H-Q, general treat- 
ments waiving those approximations did not appear un- 
til studies of Wigner functions by Haake and Reibold [5] 
who addressed the issues in low temperature and strong 
damping regimes. Exact non-Markovian master equation 
for the model was later derived by Hu, Paz and Zhang 
(HPZ) Q who employed the path-integral approach for 
initially factorizable states. The use of Wigner function 
;9|, 10 and characteristic functions [ll| presented alter- 
native means to derive the master equation. The study 
was generalized by Karrlein and Grabert Q with path- 
integral approach to cover non-factorizable initial states, 
who also pointed out that the exact Liouville operator 
for a damped harmonic oscillator is not independent of 
initial system states in general. Recently, non-Markovian 
master equations have been generalized to two-oscillator 
problems in order to study quantum disentanglement 
processes ^IUrl!- 

In parallel with these advances, there is a growing in- 
terest in coherence protection due to the need of pre- 
serving quantum information stored in a system that is 
inevitably coupled with its environment [12j . In partic- 
ular, dynamical decoupling techniques for two-level sys- 
tems based on various pul ses sequences have been de- 
veloped and optimized |13l4l7j |. For damped harmonic 
oscillator systems, parity kicks by a stream of tt pulses 
is a conceptuall y si mple strategy in realizing dynamical 
decoupling [l2|, HH . Such a strategy has been studied 
primarily in the <5-pulse regime [25l - l27j | , but little was ex- 
plored in the soft-pulse regime. In order to address the 
problem one needs to establish equations describing the 
dissipation dynamics for systems driven by time- varying 
external fields. 

The main purpose of this paper is to present an ex- 



act non-Markovian master equation for an oscillator with 
time-varying parameters in general. Such time-varying 
parameters correspond to the modulation of oscillator 
frequency and parametric interactions strengths caused 
by external devices. As a demonstrative example, we ap- 
ply our master equation to study the dynamics in 'parity 
kick' problems. Numerical results of dissipative coeffi- 
cients in the master equation are found to be drastically 
modified by the tt pulses. By solving the master equa- 
tion, we also address the effectiveness of coherence pro- 
tection, quantified by fidelity J-(t) of system state, at 
various kicking frequencies with soft pulses. 



II. THE MASTER EQUATION 

We consider an oscillator coupled linearly with an os- 
cillator bath the total Hamiltonian H — Hs(t)+HB+Hj, 
where 



Hs(t) 



Hi = J2g k (a + a>)(b k + b 



(2) 
(3) 



Here the system oscillator has a time-dependent natural 
frequency luq + f2(t) and a parametric coupling strength 
X(t). With these time-dependent parameters, we treat in 
generality the time-dependent damped harmonic oscilla- 
tor problem, which distinguishes our work from previous 
studies. Physically, the frequency shift f2(£) can be in- 
duced by off-resonance driving fields, and X(t) describes 
time-dependent parametric effects, for example in down 
conversion processes in nonlinear optics. In our model 
the bath consists of a large number of oscillators, with 
Wfe and real coupling strength g k being the frequency and 
coupling strength of the k-th mode oscillator. We also 
define the annihilation operators of the system oscillator 



2 



and the bath oscillators by a 
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and 



(^ fe m^fc ) ' Here X and P are position 



and momentum operators of system oscillator with mass 
to and frequency ujq, while q k and are position and mo- 
mentum operators of fc-th mode of bath oscillator with 
mass TOfc. 

The bath structure is characterized by spectral density 
J{lo) = X)fc5fc^( u; — I n the limit ujk becomes con- 
tinuum, we adopt the commonly used spectral density of 
the form 



J(uj) = r/Lj [ — | ( 

. OJr 



(4) 



where 77 is a dimensionless real number governing the 
strength of system-bath coupling, and u> c is the cut-off 
frequency. Necessity and justification to introduce cut- 
off frequency have been discussed in The exponent 
77 is a real number that determines the w-dependence of 
J(w) in the low frequency region, and for physical baths 
77 > 0. In literatures, < n < 1, n — I and n > 1 are 
termed as 'subohmic', 'ohmic' and 'superohmic' baths 
respectively. We will focus on ohmic bath as an example 
given in Sec. IIIII 

To derive the master equation, we first consider the 
initial total density matrix of the factorizable form: 



PtoM 



|a }(a | 



-hu k blb k /k B T 



)(1 



-hu) k /k B T\ 



) (5) 



where |«o) is a coherent state for system and the bath 
is initially in thermal equilibrium at temperature T. At 
later time t > 0, the total density matrix in position- 
space representation reads 

p tot (X",q",X',q\t) 

= J dXdq J dX"'dc{" U{X", q", t; X, q, 0) 

xp tot (X, q, X'", q'", 0)U*(X', q', t; X'\ q'", 0)(6) 



where p to t(X", q", X', q', 0) 
where q = ({qk}) and 
(X",q"|rexp \-iJ*dsH(s 



= (X",q"|p tot (0)|X',q'}, 
U(X",q",t;X',q',0) = 

\X', q'}, with T being the 



time-ordering operator. The factorized initial density 
matrix <j5j) guarantees that the Liouville operator is inde- 
pendent of initial system state, which was also observed 
in f8T-[lol]. The Guassian initial state ((5]) and the Gaus- 
sian kernel in Eq. ^ resulted from the linearity of total 
Hamiltonian allow exact integration, which makes the 
reduced density matrix ps = Trs(/3 tot ) also a Gaus- 
sian. In other words, the master equation governing 
the evolution of ps must preserve the Gaussian prop- 
erties. This requires that the master equation involves 
only some quadratic combinations of a and a) as in the 
HPZ master equations. Together with the requirements 
of conservation of probability [Tr(ps) — 0], hermiticity 



{ps = Ps) an d state-independent coefficients (71, 72, 73, 
Auj and A A in below), we have the following form of 
time-convolutionless master equation: 

ps = -i[Hs{t) + AH s {t),ps] 

-71 (t) (a f ap s + psa f a - 2ap s a f ) 

-72 (t) (aa^ps + psaa^ - 2a t p s a) 

-73 (t) (aaps + psaa - 2ap s a) 

-7|(i) (aVps + p s a)a) - 2a) Ps a)) (7) 



where 



AA(i) 2 , AA*(t) 



,t2 



+ Aw(t)a'a (8) 



are modifications to the system Hamiltonian due to 
system-bath interaction. In particular, Aw(i) is the fre- 
quency shift term, AA(t) modifies the parametric in- 
teraction, 71 (t) and 72 (t) are functions governing dissi- 
pation and amplification, and the terms with 73 (t) de- 
scribe phase-dependent decoherence typically appear in 
squeezed baths 

Our next task is to determine the time-dependent co- 
efficients. To this end, with the formal solution of bk{t), 
let us write down the Heisenberg's equation for a(t), 

a(t) = -i\*(t)a f (t)-i[uj a + n(t)}a(t) 

-iJ2dk [6 fc (0)e-^* + 6i(0) e 1 ^ 



dsK(t-s) [a(s) + a) (s)] 



(9) 



where K{t) is the memory kernel and K(t) = 
— 2i J2k 9k sm w fc r = — 2i L dcu J{u>) sinwr. 

The linearity of Eq. © leads to a general operator 
solution for a(t > 0) that can be expressed in terms of 
initial conditions, 

a(t) = G(t) a(0) + L*(t) a f (0) + F(t) (10) 

where F(t) = £ fc Uk(t)b k (0) + f£(*)&t( )) ■ Thc func- 
tions G(t), L{t), pk(t) and v k {t) can be determined by 
substituting Eq. (TT0|) into Eq. (J9j> and comparing coeffi- 
cients of initial system operators. This gives, 



G 



L = 



-i\*{t)L-i[u + tt{t)]G 

- I dsK(t-s) [G(s) + L(s)}, 
Jo 

i\(t)G + i[u} + n(t)]L 

+ f dsK{t~s) [G(s)+L(s)], 
Jo 

-iA*(t)F t - % [u + fi(t)] F 

- [ dsK{t-s) [F{s)+F\s)] 
Jo 



(11) 



(12) 



(13) 
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with initial conditions G(0) = 1, L(0) = 0, /i fe (0) = 
and ^.(0) = (see Appendix [5] for further reduction 
on F(t)). Hence for a given spectral density J(w) and 
system Hamiltonian, G(t) and L(t) can be solved and 
the operator solution a(t) can be found. 

Now we make use of the fact that the time deriva- 
tive for (a(t)), (a(t)a(t)) and (aJ(i)a(i)) obtained by the 
master equation ([7]) must agree with that obtained by the 
corresponding Heisenberg operator solution after taking 
expectation values. A direct comparison of these equa- 
tions (see Appendix iBl for details) allows us to determine 
7i> 72, 73, Aw and AA. It is easy to show in the compar- 
ison that 



AA = Auj — i("fi — 72) 



(14) 



which is not independent from other time-dependent co- 
efficients, and the remaining coefficients are: 



2iAu 



Wit) jo 



ds K(t - s) 



7i = 



272 



{[G(t)-L(t)} [L*(s) + G*(s)} 

+ [G*(t)-L*(t)][G(s)+L(s)]} (15) 

~2m + ^ (16) 

dv 1 w(ty 1 



LG* ~G*L 

W) 

L*G - GL* 



(FF) 



W(t) 



(Ft Ft) 



-273 = 



dV 1 W(t) 



(FF) 



L*G — GL* 



(F^F + FF^) 



(17) 



(18) 



where W(t) = G{t)G*(t) - L(t)L*{t) and the bath-bath 
correlation functions are 

(F(i)F(i)) = - J* ds' J ds" k t {s" - s') 

[G(s')-L*(s')} [G(s")-L*(s")](19) 

(F*(t)F(t)) = f ds' f ds"K T (s"-s') 
Jo Jo 

[G*(s') - L(s')] [G(s") - L*(s")](20) 
(F(t)F\t)) = J ds' J ds" k t {s" -s') 

[G(s')~L*(s')} [G*(s")-L(s")](2l) 

The temperature-dependent memory kernel kt(t) takes 
the form 



k t(t) = 

k 



2 cos (wfcT) e k B T — 1 



for a thermal bath at temperature T. In particular, at 
zero temperature, 



(23) 



In the special case ft(t) = X(t) = 0, we have Re (73) = 
(71 + 72)/2 and our equation can be reduced to HPZ 
master equation. 

Since the master equation ([7]) works for any initial sys- 
tem coherent state |ao), we can generalize the results 
to arbitrary initial system states by using the Glauber- 
Sudarshan P-representation. This is because for an ar- 
bitrary system state, we can formally express its density 
matrix in the diagonal form, 



p S (0)= d 2 aP(a)\a){a\ 



(24) 



By the linearity of (J7J) and the fact that all the coefficients 
are independent of initial «o, we can conclude that the 
master equation (J7J) is also valid for any initial system 
states. 



III. EXAMPLE: PARITY KICK CONTROL 

The master equation ([7]) with the coefficients given in 
Eq. (|15H18p is the main result of this paper. Such an 



equation provides a useful tool to determine the behav- 
ior of a damped harmonic oscillator subjected to external 
modulation of system parameters. To provide an explicit 
example, we employ our master equation to study the 
dynamics in 'parity kick' decoherence control problems. 
Previous studies of this subject were mostly confined to 
ideal (5-pulses or square pulses that have finite discontin- 
uous jumps [25l427| . In this section we examine parity 
kick with soft pulses and its efficiency. Specifically, we 
consider the system Hamiltonian H$(t) in Eq. (JTJ) with 
A(t) = and 



n(t) 



where 



(e,s,t) 



(e,nr + -,t) + 0(e,nr ■ 



3r 



*)} (25) 



x sin 



c 
2 

t - s 



5)] 



(26) 



(22) 



and 9(t) is the unit-step function. The driving frequency 
fl(t) corresponds to a pair of sine-squared pulses within 
one pulse period nr < t < (n + l)r, each with pulse 
duration characterized by e, and they peak at t = nr + ? 
and t = nr+^- with strengths — it and n respectively (see 
Fig. [T]). The idea of using pulses with alternating signs 
has been discussed in [2l} ■ The pair of n pulses with zero 
pulse width (e —¥ 0) can effectively reverse the direction 
of interaction. If the system and bath start uncoupled, 
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FIG. 1: Sketch of the pulses in three periods with r/2 being 
the time between successive peaks. Here e/r = 0.2 is used. 



then one can prevent the system from coupling with the 
bath at later times by applying Dirac-delta shaped pulses 
frequently. 

For ideal Dirac-delta shaped pulses (e —> 0), the pulse 
width is so short and the strength so strong as that Eqs. 
(HU and CE2]) yield G(t) -mS(t ~ t )G(t) and L(t) « 
iTr5(t — to)L(t) around t = to where the 7r-strength pulse 
peaks. Integrating, we observe that the sole effect of such 
an ideal pulse is to flip the sign of G(t) and L(t), while 
leaving G(t) and L(t) unchanged, i.e., G{t^) — > — G(£q ), 
G(to ) -> G(t+), L(to) ^ -L(t+) and L(t ) L(t+). 
This leads to the sign- flip of coefficients 71, 72, 73 and 
Aw of the above master equation. 

However, we point out that while the coefficients of 
master equation (J7J acquire negative sign immediately 
after the kick, it does not necessarily constitute an effec- 
tive scheme in coherence protection. Typically memory 
time characterizes the transient time before the coeffi- 
cients eventually settle at their long-time limits with the 
system decaying steadily. Frequent kicking leads to quick 
and repeating sign-flips that inhibits the free evolution of 
the coefficients jj(t) and their subsequent settlements, 
which also produces sawtooth-like graphs for the coeffi- 
cients. The dissipative coefficients would average to zero 
over an extended period of time and thus protect the sys- 
tem state and coherence from decaying. Specifically, for 
an effective scheme the kicking period r should be less 
than the memory time, which is of the order 1/lj c , in 
order to achieve good coherence protection. Vitali and 
Tombesi [25[ also pointed out that high-frequency bath 
oscillators with ujk 3> ujq would not be able to evolve sig- 
nificantly before interaction Hamiltonian changes sign if 
a kicking frequency much higher than cut-off frequency 
is used. Memory time, or equivalently the inverse of cut- 
off frequency, then plays an important role in fixing the 
coherence protection scheme. In light of this, the ideal 
decoupling pulse sequence would consist of pulses with 
width e — > and pulse period r <C l/u) c . More discus- 
sions about the use of 5-pulses can be found in [25M27| . 

For soft pulses, we can obtain numerically the time- 
dependent coefficients of master equation according to 
the prescription in Sec. |TTJ An example is given in Figs. 
[5] and El in which we can compare the coefficients of mas- 
ter equation under the influence of soft-pulse kicks with 



Im(A3l') 
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FIG. 2: (Colour online) Black solid line: Time-dependence 
of modification of oscillation frequency (a) and the imaginary 
part of parametric coupling (b), with 77 = 0.2, cj c /cjo = 0.2. 
e/r = 0.5 and ujqt = 2.0. Red dashed line correspond to 
cases without parity kicks for comparison. The bath is ohmic 
at zero temperature. 
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FIG. 3: (Colour online) Black solid line: Time-dependence 
of relaxation coefficients of the master equation for an ohmic 
bath at zero temperature with r\ — 0.2, oj c /^o = 0.2, e/r = 
0.5 and uqt — 2.0. Red dashed lines correspond to cases 
without parity kicks for comparison. The bath is ohmic at 
zero temperature. 



the free-evolution scenario. In our calculations, we con- 
sider an ohmic bath at zero temperature. The parameters 
77 and uj c are chosen such that non-Markovian features, 
including transient behaviour and non-exponential decay 
of coherence, can be observed for several natural oscilla- 
tion cycles. In Fig. [5]we show the modifications of system 
Hamiltonian AH$ in Eq. ([S]) due to the bath interaction. 
Without kicking (red dashed lines), Au; and A A are seen 
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to exhibit oscillatory patterns with frequency close to the 
natural frequency ujq of the system. The effect of parity 
kicks (black solid lines) seem to suppress partially both 
Aw and AA with complicated oscillations that follow the 
kicking frequency 1/r when 'parity kick' is in place. In 
Fig. [3] the time-dependence of dissipative coefficients 71, 
72 and 73 are shown. We note that perfect sign-flip is 
not present in the figures due to the use of soft pulses. 
Instead, 71 and 72 display sawtooth-like periodic oscilla- 
tions, which average to numbers close to zero in the long 
run. 
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FIG. 4: (Color online) Evolution of fidelity J-{t) subjected to 
soft pulses with cuoe — 1.0 and various kicking period r. The 
parameters are: r\ — 0.2, ui c /uio — 0.2 and an ohmic bath at 
zero temperature. 



IV. CONCLUSION 

We have derived an exact non-Markovian master equa- 
tion for a damped harmonic oscillator linearly coupled 
with an oscillator thermal bath, with a time-dependent 
natural frequency and parametric interaction strength. 
The equation enables us to determine the evolution of 
system density matrix under the influence of various 
pulses or signals. Expressions of the coefficients of mas- 
ter equation as well as the correlation functions are com- 
putable once we solve the functions G(i) and L{t). We 
emphasize that the method works because of the linear- 
ity of the whole system, which leads to the Gaussian 
propagator, so that we need only to consider a finite 
number of equations. Our approach has been applied 
to parity kick problems using soft pulses, and we observe 
that coefficients of master equations exhibit complicated 
time-dependence when the system frequency is modu- 
lated by 7r-pulses. In particular, our numerical results 
suggest that suppression of decoherence can be achieved 
by soft pulses, as long as the kicking frequency sufficiently 
higher than cut-off frequency of the bath. 
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Appendix A: derivation of G, L and F 



With the coefficients solved, we proceed to examine 
the dissipation decay of an initially excited system. As 
an illustration, we consider the initial state |^s(0)) = |1), 
where |1) is the first excited energy eigenstate of the sys- 
tem. The efficiency of decoherence protection is quan- 
tified by the fidelity F{t) = {^s$)\ps(t)Hs(G)) for 
the pure initial system state. For an ohmic bath and at 
zero temperature, we observe in Fig. 3] that free evo- 
lution (without kicks) no longer gives an exponentially 
decaying fidelity, and the initial slip is obvious. Our nu- 
merical results indicate that parity kicks with soft pulses 
can suppress the decay if r is sufficiently short. How- 
ever, we notice a transition from coherence protection 
to decoherence acceleration with increasing pulse period, 
despite the fact that the pulse periods used in Fig. [4] are 
all shorter than characteristic memory time of the bath. 
This is related to the anti-Zeno effect [28fl, and it has 
been discussed the <5-pulse limit in Ref. [27[ where the 
authors suggested that this transition, which occurs when 
t ~ l/w c , is not unique to ohmic bath. We have also per- 
formed other simulations (not shown) for soft pulses with 
various widths. For the same parameters in Fig. we 
find that increasing the width e would slightly decrease 
the fidelity. 



With the Heisenberg's equation of motion of system 
operator a(t), by combining (|9]) and (1101) we can write 
down equations (|1HI13P by comparing coefficients of op- 
erators <z(0) and a'(0). From the expectation value of 
(fTU|). we have the initial conditions G(0) = 1, L(0) = 
and (-F(O)} = 0, which would solve the differential equa- 
tions for G(t) and L(t) above given the memory kernel 
K(t). The bath operator F(t) is solved by Green's func- 
tion approach. We let 



V 



where 



f dt'T(t-t') 
Jo 



V=(F,F^Y 



B(t') 



(Al) 



(A2) 



and 



B{t) ee {B{t)-B\t)\ 



— iuikt 



MO) 



X(o) 





' 1 " 




-1 



(A3) 



The 2x2 matrix T(r) has two independent entries, which 
may be viewed as the Green functions. Differentiating 
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(TAT]) w.r.t. t and setting T(0) = -il, where I is the 2 x 2 
identity matrix, we arrive at 



[it(t - - MT(t - t') 



i dsK(t- s)T{s -t')]B(t')dt' = (A4) 



with the matrices M(t) and K(r) defined as 



M(t) = 



w + ^(r) A*(r) 
-A(r) -Lo - fi(r) 



(A5) 



equations in terms of G(t) and L(t). Expectation val- 
ues are taken with respect to the initially factorizable 
state, with the bath at thermal equilibrium. With this 
choice of initial state, system-bath correlations would be 
reduced to bath-bath correlations, e.g. 



(a(t)F(t)) = (F(t)F(t)). 



(B2) 



This follows from the fact that Tr&[pB(0)F(t)] vanishes, 
and (|B2[) is an essential condition that guarantees a state- 
independent master equation. We proceed to write down 
the equations obtained from the master equation ((7|) and 
equation (|B1|) . From the master equation ([7]), 



K(t) = K(t - s) 



1 1 

-1 -1 



(A6) 



Since the L.H.S. of equation (|A4I) must vanish for all 
B(t), we have 

it(t-t')-MT(t-t')+i f dsK(t-s)T(s-t') = (A7) 
Jo 

While one can solve the Green's functions with the initial 
condition stated above, a simple comparison with the dif- 
ferential equations for G(t) and L(t) immediately yields 
the following correspondence 



T(r) 



G(r) L*{j) 
L(r) G*(t) 



(A8) 



Thus the bath operator F(t) can be completely solved 
if coupling constants are given and G(r) and L(r) are 
known. Explicitly, 



F(t) 



i / ds[G(t-s)-L*(t-s)] 



-iui k s , j t 



b\{Q)e luJkS (A9) 



Appendix B: comparison of coefficients 

From the operator solution of a(t) (TIT)]) , taking deriva- 
tives w.r.t. t and eliminating operators at t = 0, we have 
the following operator equation 

GG* - L*L L*G — GL* + 



W(i) W(t) 
■ GG* - L*L L*G - GL* 

~r F — 



W(t) 



W(t) 



F* (Bl) 



where W(t) = G(t)G*(t) - L{t)L*(t). This operator 
equation also enables us to express second moments' 



dt 



(a) = -[at)+tx(t)}(a)-tA*(t)(a^) (B3) 



dt 



(aa) = -2[£(t) + i X (t)]{aa) - 2iA*(t)(a^a) 
-iA*(t) - 2j*(t) 



— (a f a) = -2£(t)(a?a) + iA(t)(aa) 
dt 



ik*{t){<Ja))+2 l2 (t) 



(B4) 



(B5) 



where £(t) = 71 (t) - 7 2 (t), x0) = ^0 + ^(*) + Aw(i) and 
A(t) = X(t) + AX(t). and from the Heisenberg's equation 

(EH), 



W(t) 



d , v GG* - L*Z/ , v L*G - GL* 

dt {a) = ~^v^ {a) ' 
di {aa) = 2 w(t) {aa) 

L*G - GL* , + . L*G - GL* 
+2 77777; (a^a) + 



(at) (B6) 



W(t) 



W(t) 



dt 



4{FF)-2 GG *~^ L (FF) 
dV 1 W(t) ^ 1 

L G — GL , d. ± , 

(F^F + FF^) (B7) 

aa) 



W(t) 

d . + . WYi) . t . LG* - G*L 
(a'a) = 77774( aa ) + 



W(t) 

L*G - GL 
+ - 



W(t) 



( a t a t\ + " / F t F \ 
W(t) x ' dV 1 

Wit) t LG*-G*L 

'wWV } ~ wit) {FF) 

L*G - GL 



W(t) 



(B8) 



where G and L are understood to be time-dependent 
functions. By comparing these two sets of equations we 
obtain the coefficients of master equation J7|. 



7 



[1] H. -P. Breuer and F. Petruccione, The Theory of Open 
Quantum Systems (Oxford, Great Britain, 2002). 

[2] H. J. Carmichael, Statistical Methods in Quantum Optics 
1 (Springer, Germany, 1999). 

[3] C. W. Gardiner and P. Zoller, Quantum Noise (Springer, 
Germany, 2004). 

[4] Ulrich Weiss, Quantum Dissipative Systems (World Sci- 
entific, Singapore, 1999). 

[5] Fritz Haake and Reinhard Reibold, Phys. Rev. A 32, 
2462 (1985). 

[6] Robert Karrlein and Hermann Grabert, Phys. Rev. E 55, 
153 (1997). 

[7] Chung-Hsien Chou, Ting Yu, and B. L. Hu, Phys. Rev. 

E 77, 011112 (2008). 
[8] B. L. Hu, Juan Pablo Paz, and Yuhong Zhang, Phys. 

Rev. D 45, 2843 (1992). 
[9] J. J. Halliwell and T. Yu, Phys. Rev. D 53, 2012 (1996). 
[10] G. W. Ford and R. F. O'Connell, Phys. Rev. D 64, 

105020 (2001). 
[11] Andrey Pereverzev, Phys. Rev. E 68, 026111 (2003). 
[12] Lorenza Viola, Emanuel Knill, and Seth Lloyd, Phys. 

Rev. Lett. 82, 2417 (1999). 
[13] G. S. Agarwal, M. O. Scully, and H. Walther, Phys. Rev. 

Lett. 86, 4271 (2001). 
[14] S. Pellegrin and G. Kurizki, Phys. Rev. A 71, 032328 
(2005). 



[15] P. Karbach, S. Pasini, and G. S. Uhrig, Phys. Rev. A 78, 

022315 (2008). 
[16] Gotz S. Uhrig, Phys. Rev. Lett. 98, 100504 (2007). 
[17] Gotz S. Uhrig, Phys. Rev. Lett. 102, 120502 (2009). 
[18] Leonid P. Pryadko and Gregory Quiroz, Phys. Rev. A 

80, 042317 (2009). 
[19] S. Pasini, P. Karbach, C. Raas, and G. S. Uhrig, Phys. 

Rev. A 80, 022328 (2009). 
[20] Jun-Hong An, Yc Yco, and C. H. Oh, Ann. Phys. 324, 

1737 (2009). 

[21] Jun-Hong An, Ye Yeo, Wei-Min Zhang, and C. H. Oh, J. 

Phys. A: Math. Theor. 42, 015302 (2009). 
[22] Juan Pablo Paz and Augusto J. Roncaglia, Phys. Rev. 

Lett. 100, 220401 (2008). 
[23] Ruggero Vasile, Stefano Olivares, Matteo G. A. Paris, 

and Sabrina Maniscalco, Phys. Rev. A 80, 062324 (2009). 
[24] Sabrina Maniscalco, Stefano Olivares, and Matteo G. A. 

Paris, Phys. Rev. A 75, 062119 (2007). 
[25] D. Vitali and P. Tombesi, Phys. Rev. A 59, 4178 (1999). 
[26] D. Vitali and P. Tombesi, Phys. Rev. A 65, 012305 

(2001). 

[27] K. Shiokawa and B. L. Hu, Quantum Inf. Process. 6, 55 
(2007). 

[28] A. G. Kofman and G. Kurizki, Nature 405, 546 (2000). 



